Vortex sub-lattice melting in a two-component superconductor 
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We consider the vortex matter in a three-dimensional two-component superconductor with indi- 
vidually conserved condensates with different bare phase stiffnesses in a finite magnetic field, such 
as the projected superconducting state of liquid metallic hydrogen. The ground state is a lattice 
of composite, i.e. co-centered, vortices in both order parameters. We investigate quantitatively 
two novel phase transitions when temperature is increased at fixed magnetic field, i) A "vortex 
sub-lattice melting" phase transition where vortices in the field with lowest phase stiffness ("light 
vortices") loose co-centricity with the vortices with large phase stiffness ("heavy vortices"), thus 
entering a liquid state. Remarkably, the structure factor of the light vortex sub-lattice vanishes 
continuously. This novel transition, which has no counterpart in one-component superconductors, 
is shown to be in the 3Dxy universality class. Across this transition, the lattice of heavy vortices 
remains intact, ii) A first order melting transition of the lattice of heavy vortices, with the novel 
feature that these are interacting with a background liquid of light vortices. These findings are borne 
out in large-scale Monte Carlo simulations. 

PACS numbers: 71.10.Hf, 74.10.+V, 74.90.+n,11.15.Ha 



Theories with multi-component bosonic scalar matter 
fields minimally coupled to a gauge field are of interest 
in a variety of condensed matter systems and beyond. 
This includes superconducting low temperature phases 
of light atomsi*2iMi5iii under extreme enough pressures 
to produce liquid metallic states, easy-plane quantum 
antifcrromagnetsi, as well as other multiple-component 
superconductors 2,3 . It also has applications in particle 
physics^. The projected liquid metallic states of hy- 
drogen (LMH) 1 may soon be realized in high pressure 
experiments*!*^. Moreover, LMH is abundant in the in- 
teriors of the giant planets Jupiter and Saturn, where it 
is the origin of their magnetospheresAi. At low temper- 
atures, compressed liquid hydrogen is particularly inter- 
esting since it features prominent quantum fluctuations 
which lead to the possibility of a new state of matter, 
a near ground state liquid metali. Its superconduct- 
ing counterpart involves Cooper pairs of electrons and 
protonsi, whence symmetry precludes Josephson cou- 
pling between different condensate species. Resolving 
what happens to such a system in a magnetic field is now 
a matter of some urgency, due to new and detailed first 
principles calculations predicting LMH under extreme 
pressures of order 400GPaiS. This is not far from exper- 
imentally achieved pressures of 320GPa~, where hints of 
a maximum in the melting temperature versus pressure 
are evident. Magnetic-field experiments may very likely 
be exclusive probes to provide confirmation of LMH. A 
first study of the phase diagram of the projected LMH 
in magnetic fields has been presented, unveiling a phase 
diagram with rich structure**. This raises issues of in- 
terest also in the broader domain of physics concerning 
the order and universality classes of possible phase tran- 
sitions separating phases of partially broken symmetries 
in quantum fluids. We report a study of this using a 
confluence of exact topological arguments and large-scale 



Monte-Carlo(MC) simulations. 

For general number of components N, the Ginzburg- 
Landau model is defined by the Lagrangian 
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Here, {^^'(r) | a = 1.. . N} are complex scalar fields, 
Af( Q ) is the mass of the condensate species a, and 
D = V — ieA(r). In LMH each individual condensate 

is conserved, consequently V"({^Q a ^(r)}) is only a func- 
The model is studied in the phase- 



tion of |<4 Q) (r)| 2 



only approximation 'Pq (r) 
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|^ Q) |exp[i6»( Q )(r)] where 



const. Then, the T^-term is a constant which 



may be omitted from the action^. 

For the discussions in this paper, another form of the 
action is useful. Introducing |i// Q )| 2 = \^ a) \ 2 /M^ and 

^ 2 = E a =i IV'^I 2 ! Eq. may be rewritten*! in terms 
of one charged and N — 1 neutral modes 

C = 2^2 ^J^I'VflM-e*^ +i(VxA) 2 

1 N 2 

+ 4f5 E hWhW (V(tf«-*W)) ■ (2) 

a,/3=l 

While Eq. is convenient for MC simulations, Eq. © 
has advantages for analytical considerations, since the 
neutral and the charged modes are explicitly identified. 
Moreover, Eq. @ is convenient for identifying various 
states of partially broken symmetry, emerging when an 
iV-component system is subjected to an external mag- 
netic fields. 
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We now focus on the case N = 2. In the case 
of LMH, ^/q 1 -* and \E'q 2 ' ) will denote protonic and elec- 
tronic superconducting condensates, respectively, and 
hence |^ (1) | 2 « \ip {2) \ 2 . In zero external magnetic field 
this system features a low-temperature phase-transition 
in the 3T)xy universality class at T c \ where superfluidity 
is lost, followed at higher temperatures by an inverted 
ZDxy transition at T C 2 where superconductivity and the 
Higgs mass of A (Meissner effect) is lost-. Here, we will 
consider the system in finite magnetic field at tempera- 
tures below T c2 ^. 

We define a type-a vortex as a topological defect in 
V&q^ associated with a non-trivial phase winding AO^ = 
±27r, whereas a composite vortex is a topological defect 
where type-1 and type-2 vortices coincide in space. At 




type-1 and type-2 vortices is formed, illustrated in panel 
A of Fig. At elevated temperatures, the 2-component 
system in a magnetic field will exhibit thermal excitations 
in the form of fractional-flux vortex loops similar to the 
case of zero magnetic field B = V x A = Since the 
field-induced vortices are logarithmically bound states of 
constituent (elementary) vortices, the thermal fluctua- 
tions will induce a local splitting of composite vortices in 
the form of two half-loops connected to a straight line D , 
as shown in panel B of Fig. ^ 

Consider now the processes illustrated in panel B of 
Fig. n for the case l^' 1 )] 2 <c |i/ , ^ 2 - ) | 2 , upon increasing 
the temperature beyond the low-temperature regime. We 
may view this process as a type-1 closed vortex loop su- 
perposed on a VL of (slightly) fluctuating composite vor- 
tices. An important point to notice is that a type-a vor- 
tex does not interact with a composite vortex by means 
of a neutral mode&. This follows from a topological ar- 
gument that two split branches will feature nontrivial- 
winding in the composite neutral field 6^ — 6^ 2 \ while a 
composite vortex line does not. Hence, the splitting tran- 
sition may be viewed as a type-1 vortex loop-proliferation 
in a neutral superfluid. This is illustrated in Fig. [21 
Thus, we may utilize the well-known results for the crit- 




FIG. 1: Panel A: A type-II, N = 2 superconductor at zero 
temperature in a magnetic field forms a lattice of composite 
vortices, i.e. co-centered type-1 (red) and type-2 (blue) vor- 
tices. Panel B: Low-temperature fluctuations in the compos- 
ite vortex lattice (VL). Thermal fluctuations generate closed 
loops of composite vortices and local splitting of field-induced 
composite vortices. This phase features superfluidity, as well 
as longitudinal superconductivity. Panel C: A liquid of type- 
1 vortices immersed in a background lattice of type-2 vortices. 
There is a temperature region in low magnetic field when 
type-1 vortices form a vortex liquid and the corresponding 
vortex loops are proliferated, while type-2 vortices form a VL. 
Type-1 and type-2 vortices carry only a fraction of magnetic 
flux quantum in this state. Superfluidity is lost, longitudi- 
nal superconductivity is retained. The arrow from panel B 
to panel C illustrates type-1 loop-proliferation. Panel D: A 
vortex liquid of type-1 and type-2 vortices. Superfluidity and 
longitudinal superconductivity is lost, i.e. the normal phase 5 . 
There is no type-2 loop proliferation going from panel C to 
panel D, since this is a first-order melting transition of the 
type-2 VL in the background of a liquid of linetension-less 
type-1 vortices. 

low temperatures the formation of a vortex lattice (VL) 
of non-composite vortices is forbidden because these vor- 
tices have a logarithmically divergent energy, whereas 
composite vortices have finite energy2i4. In a type-II 2- 
component superconductor, therefore, a VL of composite 



FIG. 2: Detailed illustration of the low-temperature thermal 
fluctuations in a VL of composite vortices. A local excur- 
sion of the vortex component with lowest bare phase stiffness 
(type-1 vortex) away from the composite VL may be viewed 
as a type-1 bound vortex loop superposed on the composite 
VL. The composite vortex line does not interact with a vortex 
with nontrivial winding in A7 = A(0' 1 ' — #' 2 ')-. A splitting 
transition of the composite VL, such as illustrated in going 
from panel B to panel C in Fig. may thus be viewed as a 
zero-field vortex-loop proliferation of type-1 vortices, with a 
phase transition in the 3Dxy universality class—*—*—. 

ical properties of the 3Dxy model for neutral superflu- 
ids described as a vortex-loop proliferationiSii^iM. This 
"vortex sublattice melting" phase transition is therefore 
in the SDxy universality classi2*i&ii, not a first order 
melting transition. The resulting phase is one where su- 
perfluidity is lost and longitudinal superconductivity re- 
tained in the component p illustrated in panel C of 
Fig. HJ 

Apart from the sub-lattice melting transition, thermal 
fluctuations will produce a melting transition of the type- 
2 VL at a higher temperature. It is well known that suf- 
ficiently strong thermal fluctuations drive a first-order 
melting transition of the Abrikosov lattice^ in N = 1 
superconductors. Due to the interplay with the prolif- 
erated type-1 vortices, a counterpart to this effect for 



3 



the case N = 2 when |?/' < - 1 ' ) | 2 7^ I?/''' 2 ''! 2 is more complex. 
The melting temperature Tm(B) of the type-2 Abrikosov 
lattice is suppressed with increasing magnetic field^. At 
low enough magnetic fields, upon heating the system, the 
ZDxy type-1 vortex-loop transition at T C \(B) will be en- 
countered before the melting transition of type-2 vortices 
at Tm(B). Above Tm(-B), longitudinal superconductivity 
is also lost, whence we may infer that the vortex- liquid 
mixture of liberated type-1 and type-2 vortices is the nor- 
mal metallic phased, depicted in panel D of Fig. ^ 

The above physical picture is borne out in large-scale 
MC simulations. We consider the model based on Eq. 
iJIJ for N = 2 on an L X L X L lattice with periodic 
boundary conditions for coupling constants | 2 = 0.2, 
|V> (2) | 2 = 2, ande 2 = 1/10. The ratio |^ (2) | 2 /|^ (1) | 2 = 10 
brings out one second order phase transition at T C \(B) in 
the "SOxy universality class well below the melting tem- 
perature T M {B) of the VL. In LMH |^ (2) | 2 /|^ (1) ? « 
2000, but the physical picture remains. For real esti- 
mates of T cl (B) and T M (B) in LMH, see Ref. |l. The 
Metropolis algorithm with local updating is used in com- 
bination with Ferrenberg-Swendsen reweighting 15 ;. Sys- 
tem sizes as large as L = 96 are used. The external 
magnetic field B studied is B x = B v = 0, B z = 2vr/32, 
thus there are 32 plaquettes in the (x, y)-plane per flux- 
quantum. This is imposed by splitting the gauge field 
into a static part Ao and a fluctuating part A nuct . The 
former is kept fixed to (Ag, ^(r), A§) = (0, 2ttx/, 0) 
where / = 1/32 is the magnetic filling fraction, on top of 
which the latter field is free to fluctuate. Together with 
periodic boundary conditions on A nuct , the constraint 
§ c {A-o + Afl uc t )<il = 27r/L 2 , where C is a contour en- 
closing the system in the (x,y)-plane, is ensured. It is 
imperative to fluctuate A, otherwise type-1 and type-2 
vortices do not interact^" 6 . To investigate the transi- 
tion at T c i we have performed finite size scaling (FSS) 
of the third moment of the action. The simulations are 
done by using vortices directly, similar to the procedure 
described in Ref. 0, but with a finite magnetic induction 
B z = 2tt/32. 

We compute the specific heat and the third moment 
of the action S = (3 J drC, defined as M 3 = ((S 3 ) — 
{S) 3 )/L 3 . Here, (3 is inverse temperature. The peak to 
peak value of M 3 scales with system size as 

^(l+a)/!/ and 

the width between the peaks scales as L -1 /"— . To probe 
the structural order of the vortex system we compute the 
planar structure function S'^^kx) of the local vorticity 
n ( Q )(r) = (V x [V6»( Q ) - e A]) /2tt, given by 



5 (Q) (k±) = 



1 



(fL 3 y 



(I 



(r)e 



ik 1 



I 2 ), (3) 



where r runs over dual lattice sites and k^ is perpen- 
dicular to B. This function will exhibit sharp peaks for 
the characteristic Bragg vectors K of the type-a VL and 
will feature a ring-structure in its corresponding liquid of 
type-a vortices. The signature of vortex sub-lattice melt- 
ing will be a transition from a six-fold symmetric Bragg- 
peak structure to a ring structure in SW(K) while the 



peak structure remains intact in S^(K). Furthermore, 
we compute the vortex co-centricity N co of type-1 and 
type-2 vortices, defined as iV co = iV+ — N~ Q , where 



N, 



± _ 



(r)\5 (i), . , (2), s 

v ' 1 n z (r),±n^ ' (r) 



(2) 



Ml 



(4) 



where Sij is the Kronecker-delta. The reason for consid- 
ering N co is that we then eliminate the effect of random 
overlap of vortices in the high-temperature phase T > T c \ 
due to vortex-loop proliferation, and focus on the compos- 
iteness of field-induced vortices. The quantity N co is the 
fraction of type-2 vortex segments that are co-centered 
with type-1 vortices, providing a measure of the extent 
to which vortices of type-1 and type-2 form a composite 
vortex system. Hence, it probes the splitting processes 
visualized in panel B of Fig. ^ an d in Fig. The results 
are shown in Fig. 

The specific heat has a pronounced peak at T c \ as- 
sociated with the ZDxy transition, and a broader less 
pronounced peak which is the finite field remnant of the 
zero-field inverted 3Dxy transition^. Scaling of M 3 at 
T c i shown in the inset c in Fig. [3]yiclds the critical expo- 
nents a = —0.02 ±0.05 and v = 0.67 ±0.03 in agreement 
with the SDxy universality class. A novel result, not en- 
countered in one-component superconductors, is that the 
structure function S^' (K) vanishes continuously as the 
temperature approaches T c \ from below, which is pre- 
cisely the hallmark of the decomposition transition that 
separates the two types of vortex states depicted in pan- 
els B and C in Fig. ^ A related intriguing feature is the 
vanishing in the co-centricity iV co at T c \ as a function of 
temperature, which we will discuss in detail below. The 
first-order melting transition is seen to take place at Tm, 
where S^(K) is seen to vanish discontinuously. This 
is the temperature at which the translational invariancc 
is restored through melting of the type-2 VL. The re- 
sulting translationally invariant high-temperature phase 
is depicted in panel D of Fig. In the temperature 
interval T < T c \ the system features superconductiv- 
ity and superfluidity simultaneously^, since there is long- 
range order both in the charged and the neutral vortex 
modes. In the temperature interval T c \ < T < Tm long- 
range order in the neutral mode is destroyed by loop- 
proliferation of type-1 vortices, hence superfluidity is 
losli*. However, longitudinal one-component supercon- 
ductivity is retained along the direction of the external 
magnetic field. For T > Tm superconductivity is also 
lost, hence this is the normal metallic state, which is a 
two-component vortex liquid. 

Next, we discuss these results in more detail. The most 
unusual and surprising feature is the continuous varia- 
tion of S^^K) with temperature, even at T c \ where it 
vanishes. The explanation for this is the proliferation 
of type-1 vortices (which destroys the neutral superfluid 
mode) in the background of a composite VL, which the 
type-1 vortices essentially do not see, cf. Fig. [3 As 
far as the composite neutral Bose field Q^> — is con- 
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FIG. 3: MC results for N = 2 \iP w \ 2 = 0.2, |i/> (2) | 2 = 2, 
and e = 1/vTO are presented. Panel a shows specific heat 
C v in black and co-centricity 7V co = N^, — 2V<^, in green. 
Note how the specific heat anomaly at T c \ = 0.37, associ- 
ated with the proliferation of type-1 vortices, matches the 
point at which N co drops to zero, illustrating that type-1 vor- 
tices tear themselves off type-2 vortices. The remnant of the 
zero-field anomaly in the specific heat can be seen as a hump 
at T ~ 3.6. Panel b shows the structure factors S^(K) 
in red and S 1 ' 2 ' (K) in blue for the particular Bragg vector 
K = (7r/4, — 7r/4). We see that S^^K) vanishes continuously 
at Tci, while (K) vanishes discontinuously at Tm = 2.34. 
Panel c shows the FSS plots of the M3 from which the ex- 
ponents a = —0.02 ± 0.05 and v = 0.67 ± 0.03 is extracted, 
showing that the sub- lattice melting is a critical phenomenon 
in the 3Diy universality class. Panels d, e, f, and g are gray 
shade plots of the structurefactor ^'^'(kx) for the tempera- 
tures T d = 0.35, T c = 0.4, T { = 1.66, and T g = 2.85, respec- 
tively. At Td and T e on both sides of the specific heat anomaly, 
the vortex lattice remains intact. The six fold Bragg symme- 
try is intact at Tf, but is lost at Tm to give a ring pattern at 
T g , the hallmark of a vortex liquid. 

cerned, it is precisely as if the composite VL were not 
present at all. Hence, (K) vanishes for a completely 
different reason than 5*^ 2 ^ (K) , namely due to critical fluc- 
tuations, i.e. vortex-loop proliferation in the condensate 
component with lowest bare stiffness. Such a phase tran- 



sition does not completely restore broken translational 
invariance associated with a VL, since for the type-2 
vortices quite remarkably, the VL order survives the de- 
composition transition, due to interaction between heavy 
vortices mediated by charged modes. The vanishing of 
N co is particularly interesting, and finds a natural ex- 
planation within the framework of the above discussion. 
That is, for T <C T cl , we have A^ co « 1, so the vortex 
system consists practically exclusively of composite vor- 
tices. As the temperature increases, thermal fluctuations 
induce excursions such as those illustrated in panel B of 
Fig. n and in Fig. [3 which reduces from its low- 
temperature value, reaching a minimum at T c i and then 
increase for T > T cl . Conversely, N~ a remains essen- 
tially zero until T cl , thereafter increasing monotonically. 
For temperatures above, but close to T c ±, fluctuations in 
vortices originating in Ad^ are still small, so the varia- 
tions in N co — V+ — N~ reflect thermal fluctuations in 

CO CO 

vortices originating in A6H 1 '. The increase of A^ means 
that type-1 vortex loops are thermally generated, and 
thus tend to randomly overlap more with the moderately 
fluctuating type-2 vortices. At their first order melting 
transition, type-2 vortices fluctuate only slightly. Thus, 
the vanishing of N co above T c ± reflects the increase in the 
density of thermally generated type-1 vortex loops in the 
background of a slightly fluctuating type-2 VL, cf. panel 
C in Fig. □ 

In summary, we have investigated the orders and uni- 
versality classes of thermally driven phase transitions in 
a two-component vortex system in a magnetic field. We 
find two phase transitions for the fields we consider in this 
paper (low- field regime), i) A 3Dxy phase-transition as- 
sociated with melting of the VL originating in the phases 
with lowest stiffness, which may be viewed as vortex-loop 
proliferation taking place in the background of a compos- 
ite VL. The structure function associated with the type- 

1 vortices vanishes continuously at the transition. This 
phase transition has no counterpart in a one-component 
superconductor, ii) A first-order VL melting associated 
with restoration of translational invariance of the type- 

2 vortex system. The corresponding structure function 
vanishes discontinuously at the transition, which however 
takes place in the background of a liquid of linetension- 
less type-1 vortices. This also sets the type-2 VL melt- 
ing apart from the corresponding phenomenon in one- 
component superconductors. 
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